Electron polarizability of crystalline solids in quantizing magnetic fields 
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A theory of the static electron polarizability of crystals whose energy spectrum is modified by 
quantizing magnetic fields is presented. The polarizability is strongly affected by non-dissipative 
Hall currents induced by the presence of crossed electric and magnetic fields: these can even change 
its sign. Results are illustrated in detail for a two dimensional square lattice. The polarizability and 
the Hall conductivity are respectively linked to the two topological quantum numbers entering the 
so-called Diophantine equation. These numbers could in principle be detected in actual experiments. 
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The polarizability is usually presented as an impor- 
tant property of insulators [3] . An electric field imposed 
between the plates of a capacitor penetrates into their 
interior, and charge redistribution is induced within the 
system. This effect is characterized by the polarizabil- 
ity: the ratio of induced dipole moments per unit volume 
and the local field. In metals this charge redistribution is 
non-uniform because of the charge mobility and the local 
electric field is non-uniform as well. In the particular case 
that the distribution of the atomic cores is not affected 
by the applied electric field, it is the electronic charge re- 
distribution only which determines the polarizability of 
insulators as well as metals. 

Crystalline solids have an energy spectrum composed 
of energy bands. They become insulators whenever the 
Fermi energy is located within an energy gap between 
these bands. Additional energy gaps can be opened by 
applying strong magnetic fields to two-dimensional (2D) 
as well as to three-dimensional crystals [1, H, [3| : the en- 
ergy spectrum of 2D crystals is separated in an integer 
number of subbands whenever the number of flux quanta 
contained in the unit cell corresponds to a rational num- 
ber 0] ■ The number of electrons per subband is specified 
by two topological gap numbers, one of which has been 
linked to the Hall conductance |5j. The present work 
establishes that for electrons in a "strong" periodic po- 
tential subject to a magnetic field, the remaining gap 
number is directly linked to the polarizability. We thus 
limit our attention to the static electron polarizability 
and its dependence on the electron concentration, i.e. on 
the Fermi energy fj,. This limitation allows to exclude 
dissipative processes from the consideration: no current 
flow is allowed along the applied electric field. For the 
sake of simplicity, an ensemble of spinless electrons at 
zero temperature is considered. We will also limit our 
treatment to "ideal" crystals with a rectangular unit cell 
of volume UzAq, where Qz is the lattice constant along the 
z-direction and Aq the area in the xy plane. 



In the considered geometry, the magnetic field is ap- 
plied along the z-direction, parallel with capacitor plates 
as well as with the z-th crystallographic axis. The exter- 
nal field due to the voltage drop between capacitor plates 
is applied along the y-direction. The system is assumed 
to be open along the i-direction, (i.e. periodic boundary 
conditions) allowing non-dissipative current flow. This 
condition allows the direct comparison of the polarizabil- 
ity with the topological gap numbers. Such condition 
could be realized in a Corbino geometry with cylindrical 
capacitor plates of large radius. 

We start with the case of a vanishing magnetic field. 
Generally, electrons are not equally distributed over en- 
ergy bands. The electron occupation s^f^ (/i) of the i-th 
band per volume UzAq obeys the sum rule UzAoNdi) — 

12i^f\l^) = s''°-'(^). Here N{fx) denotes the electron 
concentration (the integrated density of states). When- 
ever the Fermi energy is located within an energy gap, 
s^^\fi) = s {s integer). The wave functions of each band 
are assumed to be of Bloch-like form (i.e. extended) along 
the x-direction while along the y-direction we consider 
a Wannier-like form (i.e. localized). With this choice 
the mass-center positions of the electrons along the y- 
direction are well defined by the expectation values of the 
y-coordinate. Averaging over occupied states belonging 
to the i-th band gives an average value of the mass-center 
positions {Yi{iJ,)). An electric field along the y-direction, 
£y, leads to a redistribution of the electron charge, which 
can be characterized by the shifts of the mass-center posi- 
tions {AYi{iJ,)). These shifts are controlled by the balance 
of the electric force and of the gradient force, due to the 
background crystalline potential, which acts in opposite 
direction to bring the electrons back into their equilib- 
rium positions. Within linear response in £y, the latter 
is characterized by a force constant Ki{p,), so that: 

- sf\fi) e£y - K,{ti) {AYM = . (1) 
The resulting static electron polarizability Q!^"^(/i) is 
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given as follows 
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The above expressions are applicable even in the case of 
a non-uniform electric field along the y-direction if the 
charge redistribution does not affect the values of the 
force constants Ki. Strictly speaking, for real metallic 
systems it might be a crude approximation. 

Next, an external magnetic field splits the bands into 
subbands Because the geometry allows electrons 

to flow along the i-direction, the crossed electric £y and 
magnetic B fields give rise to a Hall current density jh- 
Consequently, an additional force - the Lorentz force - 
acts on the electron ensemble. The new balance condition 
along the y-direction becomes: 



-eN[fi)£y Jh- 
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where the index i now counts available sub-bands. The 
Hall current density has the standard form 
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where a{^) is dimensionless. Introducing s(/i) to ex- 
press the mass-center shifts {l^Yi{^)) in the form given 
by Eq. (IT|), the force balance Eq. ^ becomes: 



(5) 



where Ib = \/hc/ {eB) is the magnetic length. 

The effect of a strong magnetic field on Bloch electrons 
is more pronounced in 2D systems, which motivates our 
choice of a crystal formed by 2D planes, perpendicular 
to the magnetic field, which are separated along the z- 
direction by a^. The lattice constant az is large enough 
that electron transitions between planes are ruled out. 
We thus consider a single layer only: a 2D crystal with 
electron concentration azN{^). As was first noticed by 
Wannier and recently proved by Kellendonk Q, the 
necessary condition for the appearance of an energy gap 
in such systems reads iV(/z) = s/(azAo) -I- a /{az^-Kpg), 
where s and cr are integers, often called topological gap 
numbers, and a determines the value of the quantum 
Hall conductance [H]- The comparison of this so-called 
Diophantine equation with its general form Eq. ([5]) leads 
to the conclusion that s is related to the system polariz- 
ability. 

The relation between polarizability and Hall current is 
now illustrated with a tight-binding model on a square 
lattice, Aq = a^, assuming non-zero overlap between the 
nearest neighbor sites only. At zero magnetic field it gives 
a single energy band with cosine dispersion. The effect 



of the magnetic field is included with the Peierls substi- 
tution 0, 11]. This model was first used to obtain the 
"Hofstadter butterfly" energy spectrum For rational 
magnetic fields, satisfying: 
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(where p and i are integers), the energy spectrum is 
composed of q sub-bands well separated by energy gaps 
for which topological gap numbers are uniquely defined 

Si. 

Choosing the Landau gauge A = {—By, 0,0), the 
Hamiltonian of the system remains periodic in the x- 
direction. Zero-order wavefunctions of each row of atomic 
orbitals (j^aix — am, y — an) along the x-direction are thus 
of the following standard form: 
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(j)a{x - am,y - an) , (7) 



where is the wave number. The vector potential be- 
tween atomic sites belonging to different rows differs by 
AAj; = —BaAn, and the Peierls substitution suggests to 
shift kx entering the zero-field eigenenergies as follows: 
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so that phase factors in the overlap integrals arise be- 
tween orbitals of the neighboring rows. Assuming zero- 
order eigenfunctions in the form of a linear combina- 
tion of row-eigenfunctions, Eq. ([7]), the coefficients Cn(kx) 
have to satisfy the following Harper's equation [sl. [lo|: 
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(9) 



where Ea is the energy of the atomic orbitals (chosen to 
be zero) and AV is the overlap strength (the zero field 
band width is 8 Ay). The modulus of the c„ coefficients 
are periodic with a period q, but their amplitude differs 



by a phase j3: Cn+q = <i Cn For a given (3 G [0, 27r] 



Eq. © gives a system of q equations, and the Cn(kx,/3) 
are eigenvectors oiaqxq matrix. The eigenvalues Ep (kx) 
give the q sub-bands, each of them composed of energy 
branches determined by the phase /3, which are periodic 
in kx with period 27r/a. Because of the periodicity of the 
\cnikx, (3)\, it is natural to define Wannier-like functions, 
which are extended (localized) along the x (y) direction: 

+1 

Wn,k.Ar) = E C"+A(fc.,/3)*i°i^,,^(r) . (10) 

The average mass-center position in the y-direction reads: 
Yn,f3ikx) = / y\wn,p,kA'^fdr = an + Y^^\kx) , (11) 
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where 



(12) 



which vanishes at zero magnetic field. 
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FIG. 1: Energy spectrum Ef^ikx) (in units of tfie overlap 
strength Al/) as a function of the center of mass coordinate 
Yg {kx) (in units of the lattice constant) for the fraction 
p/q = 1/3, shown for /3 = (down triangles), P — n/2 (cir- 
cles) and l3 = n (up triangles). Empty (filled) symbols cor- 
respond to positive (negative) velocity along the x axis. The 
left and right sides show the magnetic edge states inside the 
energy gaps. 

The eigenvalue problem is mathematically equivalent 
to the one found for a weak periodic potential in a strong 
magnetic field ij], and leads to the same spectrum as a 
function of , except that the values of p and q are inter- 
changed. However, the spectrum presented as a function 
of the mass center position differs substantially from that 
of the weak periodic potential. For the ratio p/q — 1/3, 
the energy spectrum is composed of three subbands, each 
of them formed by the branches obtained for all values 
of p. The spectrum is symmetric with respect to zero 
energy because of electron-hole symmetry. This is illus- 
trated in Fig. [2 where only three values f3 = 0,7r/2 and 
TT are shown. The branches /? = 0, tt correspond to the 
subband edges, while P = n/2 characterizes the central 
branch. At the crystal edges, determined by the values 
riL and nn of the row index n, the natural condition that 
Cn^-i = Cna+1 — leads to the appearance of edge states 
formed by contributions from each of the energy branches 
p. Two types of edge states appear, non-magnetic edge 
states and magnetic ones [1, [l3| which are responsible 
for the quantum Hall effect. The latter have opposite 
velocities at opposite edges, and are shown in Fig. [T] 
for p/q — 1/3. Note that in the thermodynamic limit, 
the mass-center positions of edge states can be identified 
with the edges of the physical system. 



Next, we switch on £y, which shifts the position of 
the atomic orbitals. The force trying to return elec- 
trons back to their original positions is assumed to be 
linear in this shift, with a proportionality ("force") con- 
stant K = mof^Q, with a confinement frequency flo (mo 
is the free electron mass). Within linear response with 
respect to £y the mass-center positions of atomic or- 
bitals are shifted along the y-direction by the distance 
—e£y/{moflQ) from their equilibrium positions na. This 



shift enters the Peierls substitution, Eq. 
to an additional shift of fc^ — > fcx + A/c,. : 
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where Wc = eB/m^fC is the cyclotron frequency. The 
resulting shift of the mass-center position reads: 



^Yp(k^) = - 



dY^jk^y 



(14) 



The electric field £y also gives rise to a potential energy, 
and up to the lowest order the eigenenergies are modified 
simply by an additive contribution e£yYn^i3{kx)- As a 
result the expectation value of the velocity along the x- 
direction changes as: 



/:^vx{P,k.x) 
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(15) 



Consequently, a non-zero current density along x- 
direction (the Hall current) is induced by £y. It can be 
expressed in the form of Eq. ^ with o'(/i) given as: 
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, (16) 



where /o(-E — fi) denotes Fermi-Dirac distribution func- 
tion and N/j denotes number of branches. When the 
Fermi energy /i lies within the energy gaps, a approaches 
an integer value cr, which can differ from zero due to the 
presence of magnetic edge states. 

Using the Eq. I|14p defining the average shift of the 
mass-center positions and the identity Eq. ([S]), the elec- 
tron polarizability can be written as: 



a(/i) 



(17) 



For the tight-binding model considered here, the depen- 
dence of is fully specified by the effective topological 
number s(//). The latter is plotted as a function of the 
fiUing factor of the tight-binding band, v^, = azAoiV(^) 
in Fig. [2] for several values of the ratio p/q. The case 
p/q = corresponds to the zero magnetic field, for which 
s s*-*^^ = Vi, as follows from Eq. ([5]). For non-zero ratios 
a rich behavior is induced by the applied magnetic field. 
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FIG. 2: Dependence of the effective topological number s, on 
the band filling factor Vh for the ratio pjq = 0, 1/3 and 3/5. 

Whenever the Fermi energy is located within the energy 
gap between subbands, i.e. if I'b is given by integer mul- 
tiples of 1/g, the effective topological number s(/x) takes 
an integer value which is just equal to the topological 
gap number. As seen in Fig. [51 for p/q = 1/3, the total 
contribution of the lower and the upper sub-bands to the 
polarizability vanishes. In this particular case the applied 
electric force acting on carriers belonging to these sub- 
bands is fully compensated by the Lorentz force due to 
the induced non-dissipative edge Hall currents. For the 
case -pjq = 3/5 the magnetic field even changes the sign 
of s and consequently the sign of the polarizability. At 
half-filling of the band {vt = 1/2), the Hall current van- 
ishes and s approaches its zero-field value 1/2. The cen- 
tral symmetry of all curves around this half-filling point 
{vb = 1/2, s = 1/2) is a consequence of the electron-hole 
symmetry of the single tight-binding band. 

In this letter we have analyzed the interplay between 
the Hall current and the electron polarizability for the 
tight-binding model of a crystalline solid. This interplay 
is controlled by the compensation of the external electric 
force with the two other forces acting on the electrons: 
the gradient force, due to the periodic background po- 
tential, which is related to the polarizability, and on the 
other hand the Lorentz force, which determines the non- 
dissipative Hall current. The presence of the quantizing 
magnetic field induces a rich and complex behavior for 
the electron polarizability, which can even change sign as 
a function of the electron concentration. 

Most importantly, we provide an answer to the long- 
standing question: what measurable quantity is deter- 
mined by the topological gap number s entering the Dio- 
phantine equation? While it has been understood for 
a long time that a determines the quantum Hall effect, 
a similar interpretation of s has so far remained unclear. 
We have shown that s is directly linked to the static elec- 
tron polarizability. Contrary to the quantum Hall effect. 



the proportionality constant between the polarizability 
and the topological number s is a material dependent 
quantity rather than an universal constant, which can 
even depend on the magnetic field. The independence of 
the force constant morig on the magnetic field considered 
here is a mere consequence of the use of the Peierls sub- 
stitution which is justified when Qq ^ tOc- A shrinking 
of atomic orbitals by the magnetic field, which leads to 
a suppression of the overlap strength and thus to a de- 
crease of the overall band width, could be included by 
the substitution Op —» fig -I- lu^/A. 

The predicted effect requires that a few flux quanta 
penetrate the lattice area Aq. Even for solids with a 
lattice constant of the order of lOA, this requires exceed- 
ingly large fields 10"^ T). If, however, the magnetic 
field is tilted with respect to the high-symmetry crystal- 
lographic directions (here z), an effective area (which is 
larger than Aq) can accomodate a few flux quanta, so that 
the above condition can be reached with experimentally 
available fields This condition can also be reached 
in 2D arrays of quantum dots or antidots with a lattice 
constant ^ 100 nm. However, such (anti) dots need to 
be weakly coupled in order to reproduce a tight binding 
energy spectrum. Their size reduction (which implies a 
weak overlap and well separated atomic states) is thus 
a necessary condition for observing the predicted effect, 
and it still represents an experimental challenge. 
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